a=load('dipole_sphere_MLFMA.dat');
b=load('dipole_sphere_PO.dat');

figure(1);
plot(a(:,1),a(:,2),'b-');
hold on;
plot(a(:,1),b(:),'r+-','MarkerSize',3);
xlabel('\phi_s');
ylabel('Directivity(V/m)');
title('Vertical radiation pattern');
legend('MLFMA','PO-MM Hybrid');